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The role of quantum tunneling effect in the electron accretion current onto a negatively charged 
grain immersed in isotropic plasma is analyzed, within the quasiclassic approximation, for different 



,^3 ' plasma electron distribution functions, plasma parameters, and grain sizes. It is shown that this 

contribution can be small (negligible) for relatively large (micron-sized) dust grains in plasmas 
with electron temperatures of the order of a few eV, but becomes important for nano-sized dust 

Qh! grains (tens to hundreds nm in diameter) in cold and ultracold plasmas (electron temperatures 

1*3 

O ' tens to hundreds of Kelvin), especially in plasmas with depleted high-energy "tails" in the electron 

CZ3 

J>y energy distribution. 
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Complex plasmas - plasmas with dust particles (grains) in them (lH5| - to a large extent 
owe their complexity to the fact that the charging process of a dust grain embedded in a 
plasma is sensitive to the plasma parameters and to proximity of other grains. Understanding 
the physics of grain charging is thus important for understanding complex plasmas. 



ibrium charge of a grain immersed 
6|, |7| , in which the electron and ion 



The most commonly used model for finding the equi 
in plasma is the Orbital Motion Limited (OML) model 
currents from plasma onto the grain are found by analyzing particle orbits and determining 
whether they intersect the grain, using classical mechanics. The equilibrium grain charge is 
found from the condition that these classical currents cancel each other. However, in certain 
conditions quantum mechanical effects may become important, especially for electrons, and 
may lead to significant change of these currents, and hence of the equilibrium grain charge. 
Examples of additional currents induced by quantum effects include electron photoemission 
current (if the grain is illuminated by sufficiently energetic photons), or thermionic and/or 



field electron emission currents from the grain to the surrounding plasma 



0,y,&. 



The latter 



currents are due to quantum tunneling of electrons from a negatively charged grain through 
the grain's potential barrier into the plasma, a process analogous to that responsible for 
alpha-decay of radioactive atoms [9|, Il0| . 

The processes of spontaneous and field-assisted tunneling of thermal electrons out of neg- 
atively charged grains into plasma, and the associated emission current densities, have been 



well studied in the literature llMl5|. However, an inverse process of quantum tunneling of 



plasma electrons onto the negatively charged grain, which might, under favorable conditions, 



significantly increase the rate of electron accretion 



rom plasma onto the grain, has not re- 
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where the cross-section for electron 



ceived proper attention, with an exception of Ref. 

collisions with a spherical grain at low energies has been calculated quantum-mechanically. 

The aim of this work is thus to calculate the additional current associated with quan- 
tum tunneling of plasma electrons, that are classically forbidden to overcome the repulsive 
potential barrier, onto the negatively charged grain. We compare this additional quantum 
tunneling current with the classical electron current from plasma onto the grain, and analyze 
how this additional current affects the self-consistent equilibrium grain charge, for different 
plasma parameters and grain sizes. 

Consider a spherical grain immersed in an isotropic plasma with electrons and positive 
singly ionized ions. As the grain interacts with its environment, absorbing electrons and 



ions from the plasma and emitting electrons via processes such as photo and/or thermionic 
emission, it acquires an equilibrium net charge. Here, we assume this charge to be negative, 
which is normally the case due to higher mobility of plasma electrons compared to plasma 
ions, and relatively minor contribution of electron emission processes from plasma to the 
grain. [Note that a significant electron emission from the grain can make the total charge 
of the grain positive, in which case the electron and ion currents on the grain are found 
trivially from the classical OML theory. We will therefore not consider the case of positively 
charged grain here.] 

To determine the electron accretion current onto the grain, consider the motion of elec- 
trons in the stationary central field of the spherical negatively charged grain in isotropic 
plasma. This motion is in general described by the stationary Schrodinger's equation for 
the electron wave function ip [17J 

V 2 V + -^ [E - U{r)\ V = 0, (1) 

where E is the electron's total energy, U(r) is the electron's potential energy in the central 
field. Seeking the solution ip in form of spherical harmonics i/)(r,9,tp) = R(r)Yi m (9 , ip) , we 
obtain an equation for the radial part R(r) of the wave function, which can be easily reduced 
to the following equation for x( r ) = r -R( r ) : 



fx 

dr 2 



+ 



2m e ,^ 1(1 + 1) 



<E-U(r)) 



X = 0, (2) 



h 2 r 2 

where / is the angular momentum quantum number of the electron. This equation is equiv- 
alent to a one-dimensional Schrodinger's equation with the effective potential energy 

u ^= u ^ + ^r^= u ^ + ^ (3) 

where J 2 = h 2 l(l + 1) is the square of electron's angular momentum with respect to the 
center of the grain at r = 0. 

In what follows, we consider the quasiclassical approximation of electron's motion in 
the effective potential U e s(r), assuming that the electron's de Broglie wavelength is small 
compared to the characteristic scale of variation of U e s(r) (the criterion of validity of the 
quasiclassical approximation will be discussed below). The radial motion of an electron of 
energy E in the effective potential U e s(r) is sketched in Fig.HJ An electron with total energy 



E and absolute value of angular momentum J = hWl(l + 1), coming towards the grain from 
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FIG. 1: A sketch of radial motion of an electron incoming from infinity with energy E and angular 
momentum J in the effective central potential ([3]) of the grain of radius tq. 



U(r ).| 



classically allowed 
to hit the grain, 

EZU eff (r D ,^ 









classically prohibited 

to hit the grain, 

E < U eff (r Q , J 2 ), but can tunnel 

onto the grain with probability w t 



J2 



FIG. 2: Mapping of electron orbits (allowed and not allowed to hit the grain) on the semiplane of 
the integrals of motion E, J 2 . 

infinity, will hit the grain of radius rn if E > U e g(ro)- If E < U c s(ro), it encounters a potential 
barrier of the width n — ro, where r\ is the classical turning point of the electron defined 
from U e g(ri) = E. In classical mechanics, the electron cannot penetrate this barrier, and 
is reflected back to infinity, but in quantum mechanics the electron can tunnel through the 
barrier onto the grain with some non-zero probability Wt- In quasiclassical approximation, 
for the effective potential sketched in Fig. [H this probability can be obtained in the form 17 1 



w t 



ex Pi -- 



n 



'2ro e [U eS (r) - E]dr 



(4) 



Since the radial motion of an electron in a central field is fully defined by the two con- 
serving quantities - the electron's energy E and angular momentum J with respect to the 
center of the field, one can map all possible trajectories of the electrons, incoming from 



infinity, onto a {E, J 2 } semiplane 18| shown in Fig. |2j with the probability density of such 
trajectories defined by the distribution function of the incoming electrons expressed in terms 
of E and J 2 , / e (E, J 2 ) (here the superscript (— ) denotes the electrons with negative ra- 
dial velocities at infinity). On this semiplane, the line E = U e s(r , J 2 ) corresponds to the 
minimum energy E that an electron with angular momentum J needs to have, in order to 
be allowed by classical mechanics to hit the grain. In other words, this line separates the 
regions of parameters E, J 2 for which the electrons are allowed (above the line) or prohibited 
(below the line) by the classical mechanics to hit the grain. The corresponding "classical" 
electron current onto the grain, i.e., the current due to electrons with E > U e fi(r , J 2 ) that 



are classically allowed to hit the grain, is [18] 



41 = 4vrr 2 -^ // dEdJ 2 /H(£, J 2 ) • 1, (5) 

m e r J J E>U cS (r ,J 2 ) 

where 1 in the integral stands for the probability for an electron to hit the grain, which for 
electrons with energies E > U e s(ro, J 2 ) is equal to unity. The current on the grain due to 
electrons with E < U c s(r Q , J 2 ), which we call the "tunneling" current, is then 

41 = 4vrr 2 ^ // dEdJ 2 f$-\E, J 2 ) • w t (E, J 2 ), (6) 

m e r J J E<U cB (r ,j2) 

with the tunneling probability w t (E, J 2 ) defined by (j3j) in the quasiclassical approximation 
for electron radial motion. Note that in classical mechanics, the tunneling probability for 
electrons with E < U e g is zero [one can see this by formally taking the limit h — > in (J3J)], 
and thus I^ n = 0, as expected. With J 2 = h 2 l(l + 1) we have dJ 2 = h 2 (2l + 1), and the 
integration over J 2 in (Q turns into a sum over I: 

/£ = 4vrr 2 — Y, h2{21 + l H dE ft\E, I) ■ w t (E, I). (7) 

m e r l=Q JO 

where 

U cS (r ,l) = U(r ) + ^ l -V±^. (8) 

The tunneling probability Wt(E, J 2 ) can be evaluated analytically for a Coulomb (i.e., 
unscreened by plasma) potential of the grain, U(r) = a/r, with a = Z d e 2 , where Z d is 
the charge of the grain in electron charges. [Note that by assuming the unscreened grain 
potential, we somewhat overestimate the width of the barrier through which an electron 
with given J and E < U e g has to tunnel to reach the grain, compared to the real screened 



grain potential (e.g., the Debye-Hiickel potential). As a result, we expect to underestimate 
the tunneling probability, and thus the tunneling current obtained below for the unscreened 
grain is expected to be somewhat less than the tunneling current on a grain screened by 
plasma.] Introducing Coulomb units of mass, length, time, and energy as m e , h 2 /(am e ), 
h 3 /(a 2 m e ) and a 2 m e /h 2 , respectively, we define the corresponding dimensionless quantities 

as 

m _, m e a ~ h 2 _ / h 2 \ , . 

m = — , r = -z=-r, E = -E, n e = n e . (9) 

m e n A va e a A \m e a J 

Then, from ([7|) we have for the tunneling electron current from plasma on an unscreened 
negatively charged grain of radius r$: 

1=0 



A 2 2 °° rU eS (f ,l) 

l£=™r$> 2 (2Z + l) / dEw t (E,l)ft\E,l), (10) 

771 n , n JO 



where 



^r(ro,0 = i + ^. (11) 

fl '1 /(Z + l) ^ 1/2 
f 2f 2 



'ii 



«;*(£, = exp^-V8 / (, + ^-£) df^, (12) 



with the dimensionless classical turning point f\ of an electron with E < U e s given by 



h(E, I) = X (l+ y/l + 2E 1(1 + 1)) . (13) 



To evaluate w t (E,l) from (fT2l . we use the formulas [191 ]: 



R yR b f dx f dx 

dx = h - / — = + c 



x 2 x 2 J Xy /R J /ff' 

where R = a + bx + cx 2 , for a ^ 0, A = Aac - b 2 < 0, (14) 

V&a; + ex 2 , 2\Jbx + ex 2 f dx 

ax = he / =, tor a = 0, 15 

x 2 x J Vbx + cx 2 V ; 

with a = -E<Q,b = l,c = l(l + l)/2, and A = -(1 + 2£ /(/ + 1)) < 0. Performing the 
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FIG. 3: Dependence of tunneling probability Wt on the normalized electron energy E = E/E c 
(where E c = a 2 m e /h 2 is the energy in Coulomb units), for a fixed value of I > 0. 



integrations, we obtain for Wt(E, I) in case of unscreened grain potential: 
w t (E,l) = exp[-V8{L 2 -L 1 )], 



(16) 



r * 
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o 



1 + 2E 1(1 + 1] 



1(1 + 1) 1 + 1(1 + l)/f + y/2l(l + l)^E - 1(1 + l)/2rl - 1/ 



To 



1 + 2E 1(1 + 1) 



7T 



m _ 1 1(1 + 1) 



(17) 



(18) 



[Note that w t is a real quantity, as the imaginary parts of L\ and L 2 cancel each other.] A 
representative plot of w t (E, I) for a fixed / > is shown in Fig.[3j The probability of electron 
tunneling from the classical turning point n onto the grain increases with E until it reaches 
the value of 1 at E — U e g(fo, I) when the electron is classically allowed to hit the grain. 

The classical (J5J) and the tunneling ffTU]) electron currents make up the total electron 
current from plasma onto the grain, which charges the grain negatively. On the other hand, 
accretion of positive ions produces the ion current onto the grain which charges the grain 
positively. Since the ions are heavy and are attracted by the negatively charged grain, the 
ion current on the grain is a classical one, and is defined by the OML theory, which for the 
Maxwell-Boltzmann distribution of ions yields [2[ 

JW = V8^r 2 n t v Tl (l - e^/ty, (19) 



where r is the grain radius, rij is the ion density far from the grain, Vt% = {Ti/rrii) 1 ^ 2 is 
the ion thermal velocity, Tj is the ion temperature in units of energy, and <p s = a/ro is the 
surface potential of the grain. The balance of total electron and ion currents onto the grain 
defines the equilibrium surface potential <fr s = a/r of the grain. The corresponding equation 
of balance of currents 

41 + 41 = I {1) (20) 

is a nonlinear equation for the surface potential of the grain (and hence for the grain charge 
Zd since (p s = a/r with a = e 2 Zd), that can be solved iteratively, for a given distribution 
function / e of incoming plasma electrons far from the grain. 

Below we calculate the ratio I^i/IciL °f electron tunneling and classical currents onto 
the grain, and the ratio z c d las+tnn / Z^ of grain equilibrium charges defined from (1201) with 
and without the tunneling electron current, for several types of electron energy distributions 
and a range of plasma parameters and grain sizes. 

1. Maxwell-Boltzmann distribution of electrons 

For Maxwell-Boltzmann distribution of incoming electrons, fe (E,J 2 ) = Jm{E) = 
no e (m e /27iT e ) 3 / 2 exp(— E/T e ), where riQ e is the plasma electron density far from the grain, 
T e is the electron temperature in units of energy, we have for the classical and tunneling 
electron currents onto the grain: 

41^ = V8nrln 0e v Te exp ( —pjT J , a = e 2 Z d > 0, (21) 

Cm = ^ (27rTe)3/2 E ( 2/ + !) I dE W ^ E > l ) eM-E/T e ), (22) 



with Ues and w t defined by Eqs (TTTT) and IjTol) . respectively. 

2. Druyvesteyn distribution of electrons 

Although in theoretical calculations the Maxwellian distribution is usually assumed, the 
actual plasma electron distribution function is often significantly non- Maxwellian, such as in 
gas disharges [20(. For example, experimental measurements in inductively coupled plasma 



(ICP) [21] and capacitively coupled plasma (CCP) [22] discharges revealed that, while at 
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low pressures the electron distribution function is indeed close to Maxwellian, at higher 
pressures (p > 10 mTorr in ICP, and p > 0.5 Torr in CPP) it turns into a Druyvesteyn-like 
distribution, with depleted high-energy tail. Obviously, this should lead to a change in the 
electron classical and tunneling currents on the grain, and consequently in the equilibrium 
grain charge, compared to the case of Maxwellian electron distribution. 

For Druyvesteyn distribution of incoming electrons, /i (E,J 2 ) = /d(E) = 
n 0e A D (m e /nT e ) 3 / 2 exp[-B D (E /T e ) 2 ], where the constants A D « 0.177 and B D « 0.243 
are defined from the conditions (n e ) = riQ e and (E) = (3/2)T e (here (...) denotes the average 
over all electrons), we have for the classical and tunneling electron currents onto the grain: 



/■(e) 

'clas.D 



0.25V 8rcr n 0e v Te < -75-exp 



Ti a 
B D r Q T e 



J tun,D ~ 



0.50 



47r 2 n 0e 



1-Erf 



a 



'B 



D 



-B D 


\r„Tj 


a \ 
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r T e ) _ 


/' 



W e (27TT e )3/2 



E ( 2/ + !) 



/=o 



^eff(ro,0 



d£? Wt(E, I) exp 



*"»f 



(23) 
(24) 



where Erf is the error function, and C/ e ff and w t are defined by Eqs ffTTj) and (ITd|) . respectively. 



5. 5*iep distribution of electrons 



In complex plasmas, dust grains absorb electrons with energies above U e s(r , J 2 ), which 
depletes the high-energy "tail" of the electron distribution function [4]. For high enough 
number density of the grain component, this effect becomes significant, and the resulting 
electron distribution function can be roughly approximated by the step function 



ft\E,J 2 ) = f s (E) 



(3V2/167r)noe(me/ J B 1 




I 3 / 2 fnr P < P 

for E > E m »*. 



(25) 



Here, E max = e\4> s \ + 5E, where e\4> s \ is the minimum energy at which an electron with / = 
(zero angular momentum with respect to the grain) hits the grain, and 5E accounts for the 
fact that most electrons have / > 0. For simplicity, SE can be approximated by a constant, 
defined, e.g., from experimental data for a particular complex plasma system. Here for 
definitiveness we will assume a reasonable value SE = 0.5T e . For the step distribution of 
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FIG. 4: Comparison of i]i = itun/-vus as a f unc ti° n of the grain radius ro, for different electron 
distribution functions: Maxwellian (solid line), Druyvesteyn (dashed line), and the step distribution 
with 5E = 0.5T e (dash-dotted line). The currents are calculated for grains with corresponding 
equilibrium charges, immersed in argon plasma with T e = 0.1 eV and T e /Ti = 10 2 . 



incoming electrons, we have for the classical and tunneling electron currents onto the grain: 



r(e) 

'clas.S 



3tt 



--r 2 Q n 0e 



5E 2 



2^Mr e , ^\ a / ro + 5E f 2 ' 



( e ) 3vr / 2 n 0e a' 

tun ' s : : 4 V m e p 3 /2 



00 fUeS(f ,l) 

V (2/ + 1) / d^ w 4 (E, l)a H (E max - E), 

Jo 



(26) 
(27) 



^rnax ; =0 

where <r^ is the Heaviside step function, and t/ eff and w t are denned by Eqs (fTTj) and ([TBI) , 
respectively. 

The importance of the electron tunneling current onto the grain can be characterized by 
the ratio of the tunneling and classical electron currents r\i = Hun/^dL- The tunneling elec- 
tron accretion current increases the total electron current onto the grain, thus offsetting the 
total electron-ion current balance, and changing the equilibrium grain charge. This change 
is characterized by the ratio r/z = ^ las+tun /^ las of the grain charges defined from (1201) with 
and without accounting for the electron tunneling current. The ratios t]i and r\z depend on 
several factors: the electron distribution function, electron and ion temperatures, the grain 
size, and the type of gas used in the plasma discharge. These dependencies are illustrated in 
Fig. H] and in Tables |T] and [TTJ Generally, since r\i is proportional to the ratio of populations 
of electrons with E > U c s(r , J 2 ) and E < U c s(r , J 2 ), both r\i and r\z increase for electron 
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TABLE I: Ratio ^tun/^cUs °^ tunneling and classical electron accretion currents onto the grain, 
for different electron distribution functions, electron temperatures, and grain sizes. The currents 
are calculated for grains with corresponding equilibrium charges, defined by Eq. ()20p . immersed in 
argon plasma with T e /Ti = 10 2 . 
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TABLE II: Ratio zf as+tan /Z^ aa of equilibrium grain charges defined by the current balance (f20|) 
with and without the electron tunneling current, for different electron distribution functions, elec- 
tron temperatures, and grain sizes. The charges are calculated for grains in argon plasma with 
T e /Ti = 10 2 . 

distributions with depleted tails, such as the Druyvesteyn or the step distribution, as seen 
in Fig. HI Due to the same reason, for a given type of electron distribution function, rji and 
r)z increase at lower electron temperatures T e (for a fixed ion temperature Tj), and increase 
for smaller grain sizes r , as seen from Fig. H] and Tables [U and [TTJ [Remarkably, since the 
ratio of populations of electrons with E > U c s(r , J 2 ) and E < cV eff (r , J 2 ) increases with 
the slope of the E = U c s line in Fig. [2] (i.e., decreases with the grain size r ), the dependence 
of rjj on the grain size r can be well approximated by a universal power law, r]j oc Tq ' 7 , 
where the coefficient of proportionality only depends on the plasma parameters (i.e., elec- 
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tron distribution function, electron and ion temperatures, gas atomic mass), and is (almost) 
independent of ro.] 

We should note that actually the electron distribution function almost never has a clear- 
cut form such as Maxwellian or Druyvesteyn or step function, but instead can be a superpo- 
sition of several distributions of different type with different effective temperatures. Hence, 
the results of calculations of rji and rjz shown in Fig. H] and in Tables HHIT1 should be perceived 
as illustrating the trends. But if the electron distribution function can be measured, the 
corresponding electron current onto the grain can be calculated using Eqs (jSJ) and (J7J) with 
this distribution. 

The criterion of validity of the quasiclassical approximation, in which the expression (j3j) 
for the tunneling probability is obtained, reduces to the requirement that the electron's de 



Broglie wavelength X ~ h/ y/2m e \E\ is small compared to the size a/l-E 1 ] of the region near 

171 | . Because C/ e g is minimum for 



the grain where the electron energy E is of the order of U c g 

the electrons with / = 0, this criterion is the strongest for such electrons, i.e., if it is satisfied 

for electrons with I = 0, it is automatically satisfied for electrons with I > 0. For I = 0, this 



criterion reduces to hv/a <C 1, where v ~ \J ' E /m e is the classical electron velocity. With 



v ~ Vxe — yT e /m e and a ~ r zT e , where z = e\<f> s \/T e is the dimensionless grain charge, 
z ~ 1, we have the following requirement for the grain size r for which the quasiclassical 
approximation of electron motion is applicable: 

r » - 7 t=. (28) 

For T e ~ 1 eV this gives r ^> 1CT 8 cm, and for T e ~ 0.01 eV this gives r ^> 3 ■ 10~ 7 cm, 
which is well satisfied for the grain sizes used in Fig. H] and Tables U and [TT1 We note 
that for very small grains, or for very low electron temperatures, when the quasiclassical 
approximation becomes invalid, the Schrodinger's equation fl2]) for the radial wave function 
of an electron has to be solved exactly, in order to define the electron current from plasma 
to the grain. This is however hardly necessary, as the criterion ( 128]) is well satisfied for a 
wide range of plasma electron temperatures and grain sizes. 

As seen from Fig. H] and Tables [J and [Til the effect of electron tunneling from plasma 
onto the grain is most pronounced and significant for small grains (ro ~ tens to hundreds 
of nm) in plasmas with low electron temperatures (T e ~ tens to hundreds of K), especially 
for electron distributions with depleted high-energy "tails". Note that we underestimated 
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the tunneling current by assuming the grain field to be unscreened. In reality the grain is 
shielded by plasma and its potential decays faster than the Coulomb potential, narrowing 
the width of the potential barrier around the grain, and thus further increasing the electron 
tunneling probability and the electron accretion current onto the grain, especially in low 
temperature plasmas with small Debye lengths. We therefore expect the considered effect of 
quantum tunneling of plasma electrons onto the grain to be important in plasmas with nano 



0,3, 



and submicron size dust grains |l|, |3| , in ultracold plasmas (where the electron 



temperature 



can be as low as 30 K [23j), and in dark molecular clouds in astrophysics [ljj |24 |. 
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